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■ Abstract 

t-H ! 

We present a review of atmospheric and underground muon flux mea- 

| surements. The relevance of these data for the atmospheric neutrino flux 

computation is emphasized. Possible sources of systematic errors in the mea- 

^ ■ surements are discussed, focusing on the sea level muon data. Underground 

O ■ muon data are also reported. 

O ■ 

o : 

* ! 1 Introduction 

V ■ 

At sea level, together with neutrinos, muons are the most abundant parti- 
^ . cles originated by the interactions of primary cosmic rays at the top of the 

>- . atmosphere. Due to their relative stability and small cross sections, these 

particle are able to arrive deep underground and/or deep underwater. As a 
consequence, their study covers many aspects of cosmic ray physics. 

Recently, atmospheric muon flux measurements received attention the 
context of the atmospheric neutrino anomaly [I], Because of the close 
relation between muon and neutrino production, it follows that the evalu- 
ation of the atmospheric muon flux can provide an important cross check 
on the atmospheric neutrino flux. Moreover, measurements of muon flux at 
all geomagnetic latitudes are crucial for the normalizazion of the calculated 
neutrino flux. 

Measurements performed at different altitudes (sea level, at mountain 
level or in balloon born experiments) offer various advantages. First of all, 
a different set of measurements at different altitude values can provide in- 
formations about the longitudinal development of the muon component in 
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cosmic ray showers. Moreover, the interpretation of data collected at the top 
of the atmosphere are not affected by the uncertainties inherent in particle 
production and propagation, since the muons (and the corresponding neu- 
trinos) generated in the first stages of the cascade. Finally, the knowledge 
of high altitude muon data is crucial for sea level sub-GeV neutrinos: the 
corresponding sub-GeV muons originated in the same decay processes can- 
not reach the sea level, considering that the average muon energy loss in the 
atmosphere is of the order of 2 GeV; we are thus forced to take data at high 
altitudes. 

On the other hand, measurements performed at ground level offer the 
advantage of a high stability, large collecting factor and a long exposure time 
due to the relatively favourable experimental conditions. They however suf- 
fer of an intrinsic difficulty in interpretation, since the muons that arrive 
at sea level are the last stage of a multi-step cascade process. This is true, 
in particular, for the measurements at high zenith angles, near the horizon, 
where the intermediate and high energy regions of the spectrum can be anal- 
ysed (pfj,— 10-100 GeV/c). Nevertheless, for this reason, sea level data offer 
the possibility to perform a robust check of the reliability of existing Monte 
Carlo codes. 

Finally, underground measurements offer the possibility to extend the 
energy range of the muon spectrum beyond 1 TeV. Such measurements are 
of an indirect type, but their link with the direct low-energy observations 
gives the possibility to complete the picture of muon spectra measurements 
and to cross-check the validity of the global set of data. 

Most of the experiments devoted to the measurement of the muon mo- 
mentum spectra and intensity have been carried in the '70s. The problem is 
that the results are often in disagreements with one another; the discrepan- 
cies are significantly larger than the experimental errors reported. Recently 
new instruments, mainly designed for balloon experiments, have been devel- 
oped; they are able to give detailed information on the muon flux at different 
altitudes in the atmosphere 0. Also new measurements deep underground 
j|, §, §|, 0, |J or by EAS arrays || have added new information at very high 
energies. 

In this paper we summarize the observations of the muon flux at sea level 
and deep underground and discuss some of the systematics connected with 
such measurements. For more complete discussion one can refer to the recent 
papers 



10, 111 and to books 
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2 Atmospheric muon production and propa- 
gation 

Secondary muons are mainly produced in the decays of secondary mesons, 
mostly ^ and K ± . The most important decay channels, and their respective 
decay probabilities, are: 

a) tt ± -> / u ± z/ M ~ 100% 

b) K± -> v^Vp ~ 63.5% 

in which the produced muons take on the average 79% and 52% of the energy 
of the 7r ± and K ± , respectively. 

The contribution of K decays to muon production is a function of the 
energy and ranges from ~ 5% at low energies to an asymptotic value of ~ 
27% for E > 1 TeV. At very high energies a small contribution arises from 
charmed particles. The analytical form of the muon production spectrum 
at a given height in the atmosphere can be derived by folding the two-body 
decay kinematics of the parent mesons with their production spectrum. The 
latter is generally expressed in terms of the so called "spectrum weighted" 
moments 

Z^ = f\i^dx (1) 

where dN p7r ±/dx is the pion production spectrum (x = E n /E p and 7 is the 
differential primary spectral index). A similar expression can be obtained for 
kaons. In general, the development of the meson and muon components in 
the atmosphere depends on the energy range we are considering. The com- 
petition between interaction and decay of the particles plays a crucial role 
and the relative importance of the two processes depends on the energy. We 
can distinguish three different energy regions in the muon spectrum: 

a) ^> e nj K ) where e n = 115 GeV and ex = 850 GeV are the crit- 
ical energy beyond which meson reinteractions cannot be neglected. This 
is the typical muon energy range studied by underground detectors or by 
ground based experiments looking at high inclined directions. In this case, 
the meson production spectrum have the same power law dependence of the 
primary cosmic rays, but the rate of their decay has an extra E^ 1 depen- 
dence with respect to the primary and meson spectrum (a consequence of 
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the Lorentz time dilatation). The muon (and hence neutrino) flux takes the 
form: dN/dE^ = E~^ +1 \ and a zenith dependence dN/dcos6 oc (cosO)^ 1 . 
It should be noted that, in this energy region, the enhancement of the 
contribution to the secondary lepton production is particularly important in 
the neutrino flux calculation, as a consequence of the two-body decay kine- 
matics |14]]. This last remark does not hold for muons, for which the limited 



knowledge of meson production in this energy range is not so crucial as for 
neutrinos. 

b) < < e nt K , where e M ~ 1 GeV. In this energy range, almost all 
the mesons decay, and the muon flux has a power law dependence with the 
same spectral index of the parent mesons (and hence of the primary cosmic 
ray, in the assumption of complete Feynman scaling validity) and is almost 
independent on the zenith angle. A compact form which expresses the low 



and high energy regions is [12]: 
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c) E^ < e M . In this case, muon decays and the energy losses in the 
atmosphere cannot be neglected. Moreover, geomagnetic latitude and solar 
modulation now play an important role being the primary cosmic ray energy 
E p < 20 GeV. 

We stress again the relevance of muon flux measurements for the knowl- 
edge of the neutrino flux. In principle, sea level neutrino flux computation 
can be derived directly from muon flux measurements high in the atmosphere 
(X < 37 g/cm 2 ) Jl5|]. This approach gives good results, but only a complete 
Monte Carlo simulation can take into account second order effects. The main 
ingredients in Monte Carlo calculations (and the main sources of systematics) 
of atmospheric lepton production are the input primary cosmic ray spectrum 
and a detailed description of secondary multiparticle production in the at- 
mosphere. The primary cosmic ray composition plays an important role only 
in the very high energy range; the composition is dominated by protons and 
a particles at energies below 100 GeV. 

Among various Monte Carlo codes now available, we recall |l| and [16 



which are the ones used to interpret the atmospheric neutrino anomaly, and 



the new code based on the FLUKA interaction model which takes into 
account 3-dimentional effects of secondary propagation in the atmosphere. 
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The comparison between the Monte Carlo evaluation of the muon flux at 
different altitudes and at sea level with the existing muon flux measurements 
constitutes one of the most powerful benchmark to assess the validity of the 
simulations. 



3 Atmospheric muon flux measurements 

Measurements of the absolute intensity, energy spectrum and positive-to- 
negative ratio of muons have been carried out many times in the past. Most 
of these observations were made at sea level and few at different mountain 
altitudes with counter telescopes separated by absorbers (Pb, Fe) and mag- 
netic spectrometers. More recently, with the development of superconducting 
magnet, it has been possible to operate spectrometers also on board of bal- 



loons |18], [19], |2(], [21], |22| which led to accurate measurements at different 
levels in the atmosphere. 

Here we will consider mainly ground-based measurements and those made 
with detectors on balloons near the ground level or very close to it. The 
relevant quantities that can be directly measured and will be discussed here, 
are: 

- absolute muon intensity 

- muon momentum spectrum 

- charge ratio 

3.1 Absolute intensity measurements 

The vertical muon intensity at sea level is a quantity which varies with the 
geomagnetic latitude, altitude, solar activity and atmospheric conditions. 

The geomagnetic field tends to prevent low energy cosmic rays from pen- 
etrating through the magnetosphere down to the Earth's atmosphere. At 
any point on the Earth one can define a threshold or cut-off rigidity, Pc, for 
cosmic rays arriving at a particular zenith and azimuth angle [T2^ . 

Primary nuclei having lower rigidity are excluded by the action of the 
geomagnetic field and do not contribute to production of secondaries in the 
atmosphere. The cut-off values range from less than 1 GV near the geo- 



magnetic poles to about 16 GV for vertical particles near the equator [24 



It results that geomagnetic effects are important for sea level muons up to 
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about ~ 5 GeV (Fig. |1]). The effect is larger at higher altitudes; Conversi p5| 
found that the vertical flux of muons with momentum around 0.33 GeV/c at 
latitude 60 deg was 1.8 times higher with respect to the flux at the equator. 
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Figure 1: The latitude effect on the integral muon intensity at sea level. It is 
shown the ratio between the intensity measured at Kiel Ik (Pc = 2.3 GV), 
and the intensity measured near the equator I eq (Pc = lJ^.l GV), with the 
same instrument 26 . 

Moreover, as cosmic ray primaries are predominantly positively charged 
particles, the flux and spectra in the East and West directions differ up to 
energies of about 100 GeV; the intensity from the West is stronger than that 
from the East. This effect increases with altitude. 

In addition, the primary cosmic ray spectrum at the top of the atmosphere 
changes with the 11 year solar cycle as the configuration of the Interplanetary 
Magnetic Field (IMF) varies. It results that the cosmic ray flux is significantly 
"modulated" up to energies of about 20 GeV (Fig. g). 
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Figure 2: Approximate solutions (the so called force-field solution) 27,28 that 
fit the observations of the primary proton spectrum at the top of the atmo- 
sphere at different years. Also shown is the assumed Interstellar Spectrum. 
Notice the extension of the solar modulation effect. 



In order to estimate how these changes in the primary spectrum influence 
the counting rate of a muon detector it is necessary to know the " differential 
response curve" |I2"|] . Its shape varies significantly with the depth of obser- 
vations, see Fig. [j| Their detailed calculations depend on the properties of 
nuclear cascades in the atmosphere; more precise descriptions can be found 



in [RO, pTJ At the standard momentum of 1 GeV/c and at high latitudes 



the modulation is 7% and 4.5% for the differential and the integral fluxes, 



respectively |32| . 

So in making a comparison of muon observations at low energies (less 
than 20 GeV) it is very important to know the year and the location the 
measurements were made. Figure El shows the neutron monitor counting rate 
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recorded by a middle latitude station since 1953. No continuous recording 
of the same kind exists for muon monitors. By the comparison of the peak 
to peak variations during the interval 1965-1972 one can estimate that the 
total muon flux changes are usually a facto 3 to 5 smaller than the observed 
neutron flux variations P2L K31 

Finally, changes in pressure and, particularly, temperature above the in- 
strument up to the point of muon production by pions and kaons, produce 
variations of different amplitude in different energy range. The most con- 



spicuous for muons at higher energies is the seasonal variation [[3J, |35|, [36 



for which the results reported in the following have not been corrected for. 

Classical definition of the hard component [37] is related to penetration 
characteristics, it is to say the capability of crossing 167 g/cm 2 , equivalent to 
roughly 15 cm of Pb. As a matter of fact this component is made of of muon 
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Figure 4: Observed time variation in the monthly counting rate of the Climax 
Neutron Monitor (Pc = 3.03 GV) normalized to 1965. The 11 year variation 
associated with solar sunspot cycle is easily recognizable. Pronounced min- 
ima occur during sunspot maxima. A complete record of the muon intensity 
variations during the same period is not available. 



with momenta p^ > 0.32 GeV/c and less than 1% are protons, neutrons, 
electrons and pions. 

Let us distinguish between vertical and horizontal integral muon fluxes. 
These latter are made in order to extend the range of the former beyond 
several tens of GeV/c, but usually they do not give absolute values of the 
vertical muon intensity. For this reason we will not discuss them here. 

The first measurement of the integral vertical intensity was made by 
Greisen |38[ at latitude 50° and altitude 259 m a.s.l. (corresponding to 
1007 g/cm 2 ) who found the value: 0.83 x 10~ 2 ± 1% cmT 2 s _1 sr _1 . Rossi 
[ j37j] noticed that this value needed to be corrected in order to account for 
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showering and scattering of particles inside the apparatus. Successive mea- 
surements led to higher values. The observations made at different latitudes 
and during different years are presented in Fig. |5|. Most measurements were 

and only the few at low lati- 



made at high latitudes [39, HO, |41| 
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tudes [H5| , |4q| ; were corrected for the geomagnetic effect. No corrections have 
been made for the solar modulation effects; the measurements are essentially 
grouped in the period 1967-1977, with one in 1998 [44]. The agreement be- 



tween the measurements is fairly good (all the data within 10%) and one has 
to take into account that the largest contribution to the deviations are the 
systematic errors due to incorrect knowledge of the acceptance, efficiency of 
the counters and correction for the multiple scattering. 
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Figure 5: Integral momentum spectrum of low energy muons at sea level. 
Data are taken from Kiel 26 > 32 ; Flint et al. 39 , Karmakar et al. 45 , Baschiera 
et al. Al , De et al. 46 , Ng et al. 40 ; Ashton et al. 42 , Barbuti et al. 
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For energies En > 1 TeV direct measurments of the muon flux were made 
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at highly inclined directions using large magnetic spectrometers fl7\, ^8], 
large emulsion chambers [50 and EAS arrays 0. For the discussion of the 
results we refer to the recent paper ITT . 



3.2 Momentum spectra 

These spectra have been measured many times for moments up to ~ 100 
TeV/c. Magnetic spectrometers are mainly used at low and intermediate 
energies, while observations at high energies are made close to the horizon- 
tal directions at ground level or deep underground. The latter are indirect 
measurements, since the ground level spectra have to be extracted from un- 
derground data. We will consider here only ground level and underground 
observations. 



3.2.1 Ground level measurements 

Direct measurements of momentum spectra for < 1 TeV/c are important 
for the comparison of nuclear cascade models with available data. Further- 
more by extending the model results to higher energies one can hope to be 
able to evaluate prompt muon production and/or charm production. In the 
momentum region 10 GeV/c - 1 TeV/c: a) the production spectrum of the 
charged pions cannot be represented by a power law but has a maximum at 
an energy that depends on both the altitude and the latitude b) the energy 
loss and the decay of muons must be properly considered. In order to join 
low to high momentum spectra it is important to have single experiments 
that cover the widest energy range. 

To better see the differences between the sea-level spectra we plot in Fig. 
|^ the percentage deviations of the data from the best fit spectrum obtained 



by |5T |. Notice that even if individual errors are small (however increasing 
with momentum due to decreasing number of detectable particles and to 
the maximum detectable momentum), deviations up to ± 20% are observed 
probably because of systematic effects. 



3.2.2 Underground measurements 

From underground muon intensity measurements, informations about sea- 
level muon spectra can be obtained using different procedures (see for ex- 
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Figure 6: Relative deviation of the differential muon spectra with respect to 
the Kiel fit 51 . The data are taken from Allkofer et al. 51 , MASS 44 , Ayre et 
al. 52 , Bateman et al. 53 , Barber et al. 
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ample || ||). Here and in the following, we assume a standard procedure 
applied form large area underground experiment ]8|, [7|. The vertical muon 
intensity, for a given direction 6, <fi and a corresponding rock slant depth h 
can be expressed as: 

tffr = (*r) v AO f (3) 
p VAT/ / cos Oj 

where AT is the total livetime of the experiment, iVj is the number of 
detected events with multiplicity in the angular bin AQj, Aj and ej are, 
respectively, the geometrical and intrinsic acceptance of the detector. The 
relation between the measured I^(h) and the sea- level muon spectrum can 
be expressed as: 
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where P(E, h) is the muon survival probability function determined via 
Monte Carlo. Assuming for the sea-level muon spectrum an expression of the 
form (U), leaving as free parameters the muon spectral index and a normal- 
ization constant, it is possible it is possible to unfold sea level muon spectrum 
from the measured absolute muon intensity. 

In Fig. |7] are reported the results of the fit of MACRO data || together 
with LVD 0, MSU and Baksan || data. Data are presented multiplied by 
factor p 3 to better observe the variation of the spectrum in the whole energy 
region and to strengthen a possible flattening in the tail of the spectrum due 
to charm production. The statistics is still too poor to allow any definite 
assessment on the existence of this effect at energies > 10 TeV/c. 

In indirect measurements, accurate estimates of the systematic errors are 
needed. The main sources of systematics in (|]) are the knowledge of the 
rock density overburden and the treatment of hard processes in the energy 
loss of muons in the rock. In the MACRO fit, for example, their overall 
contribution has been estimated to be ~ 5% and 3% in the determination of 
the normalization constant and muon spectral index respectively. 



3.3 Charge ratio 

In the primary cosmic rays there is an excess of positively charged particles 
(protons) with respect to the total number of nucleons. This excess is trans- 
mitted via nuclear interactions to pions and further to muons. By assuming 
that the primary composition is constant in the energy range considered, 
this ratio will remain constant with the exception of high energies, where the 
contribution from kaons starts to become sizeable. The muon charge ratio is 
expected to increase also with zenith angle as the depth is increasing and like- 
wise the energy of the primaries that produce muons of a given momentum 
at ground. This quantity is important to study nucleon-nucleon interactions, 
composition and kaon contribution. Magnetic spectrographs are used for 
determining this ratio. Because of systematic effects in the momentum mea- 
surement the values are usually much spread out. Moreover limited statistics 
at high energy makes it difficult to appreciate the energy dependence. We 



report in Fig. |8] only the recent data from Mass at lower momenta, 
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Figure 7: Vertical differential momentum of muons at sea level. Direct data 
extend up to 1 TeV and are taken from [allkofer, depascale, barber, bateman, 
ayre]. Indirect data obtained from underground observations for E>1 TeV/c 
are taken from Baksan 5 , LVD 7 , MSU 6 , and MACRO fit 8 , with the two 
parallel lines indicating the range allowed by the errors of the fit parameters. 



and the two compilations made by |j55 |. It is clear more measurements with 
longer exposures are still needed. 



4 Conclusions 



Since atmospheric muons and neutrinos are generated in the same processes, 
the accuracy of the neutrino flux calculation can be improved by forcing the 
poorly known input parameters of the cascade model to fit the data on the 
muon flux. 
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Figure 8: i?atzo at sea level. The data have been taken from MASS 44 , 

Rastin 55 . 



However, the data are still not sufficient for this purpose, since several 
sea level measurements of the vertical muon flux are in poor agreement with 
one another, even though each experiment has typically very good statistics. 

Disagreement between the results of different experiments are present 
even if the quoted errors are relatively small in the majority of the exper- 
iments. It indicates the existence of significant systematic errors in some 
experiments by as much as 30-35% at momenta from 10 to 1000 GeV/c. 
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